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In  order  to  validate  existing  techniques  for  computer 
simulation  of  the  welding  process,  an  experiment  is  designed 
and  carried  out  in  which  seam  and  girth  welds  are  produced 
on  a  high  strength  steel  cylinder.   Fabrication  of  the  HY-130 
cylinder,  selection  of  welding  parameters  for  the  seam  and 
girth  welds,  and  the  selection  and  arrangement  of  instrumenta- 
tion are  described.   The  theoretical  background  of  finite 
element  computer  simulation  of  heat  flow  in  weldments  as  well 
as  techniques  for  axisymmetric  finite  element  analysis  of 
thermal  stress  and  metal  movement  during  welding  are  dis- 
cussed.  Finite  element  representation  of  weldment  geometry 
for  heat  flow  and  stress  analysis  is  presented.   Metallurgical 
properties  of  the  HY-130  cylinder  seam  weld  are  presented 
in  the  form  of  photo-micrographs  and  electron  scan  micro- 
graphs . 

Temperature,  strain,  and  radial  distortion  are  measured 
at  various  locations  on  the  HY-130  cylinder  and  the  measured 
values  are  compared  in  tabular  and  graphical  form  with  those 
values  predicted  by  use  of  the  heat  flow  and  stress  analysis 
programs.   Excellent  agreement  is  achieved  between  measured 
and  analytically  predicted  temperature  distributions,  but 
only  qualitative  agreement  is  shown  for  strains  and  radial 
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CHAPTER  I 

I. A   INTRODUCTION: 

Welding  has  been  used  as  the  primary  means  of  fabricating 
ocean  structures  for  the  past  forty  years.   Design  tools  for 
the  fabrication  of  welded  structures  have  evolved  primarily 
as  the  result  of  experimentation.   These  tools  are  usually 
expressed  as  empirical  relationships  among  welding  parameters, 
some  of  which  include  structure  and  weld  joint  geometry, 
physical  properties  of  the  structural  materials  and  welding 
consummables ,  as  well  as  methods  and  rates  of  energy  trans- 
fer during  the  joining  process.   In  most  cases  the  use  of 
these  empirical  relationships  has  been  adequate  to  ensure  the 
achievement  of  good  quality  weld  joints  in  a  production  en- 
vironment. 

As  design  tools  become  more  sophisticated,  and  the  in- 
terrelationships among  physical  properties  of  welding  and 
structural  materials  are  better  understood,  design  practices 
relating  to  allowable  stresses  in  welded  structures  converge 
toward  material  limits.   Development  and  refinement  of  tech- 
niques for  detecting  and  classifying  flaws  in  weldments  is 
another  design  consideration  which  tends  to  increase  confi- 
dence in  welded  structures  subjected  to  more  demanding  re- 
quirements, allowing  design  engineers  to  increase  allowable 
stresses  or  to  reduce  factors  of  safety  or  uncertainty. 
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An  area  of  welding  technology  which  requires  further 
definition  and  refinement  is  the  occurrence  of  residual 
stresses  and  structural  distortion  during  and  as  the  result 
of  the  welding  process.   Extreme  non-uniform  thermal  cycles 
of  weld  metal  and  of  structural  metal  in  areas  near  the 
weldment  produce  unevenly  distributed  thermal  stresses 
which  cause  localized  plastic  deformations,  with  resultant 
residual  stresses  and  structural  distortion.   These  residual 
stresses  can  significantly  contribute  to  structural  failure, 
particularly  when  the  failure  mechanism  is  dependent  upon 
low  applied  stress,  such  as  brittle  fracture  or  failure  re- 
sulting from  stress  corrosion  cracking.   Distortions  caused 
by  welding  can  also  cause  structural  failure  by  reducing  a 
structure's  ability  to  withstand  buckling  loads. 

The  ability  to  accurately  predict  the  magnitude  and 
type  of  residual  stresses  and  distortion  resulting  from  the 
welding  of  ocean  structures  would  greatly  enhance  a  design 
engineer's  ability  to  estimate  total  stresses  which  materials 
must  withstand  during  loading.   More  accurate  estimation  of 
total  stresses  would  allow  for  more  effective  use  of  materials 
and  the  increased  predictability  of  structural  failure  would 
ultimately  result  in  the  production  of  more  reliable  ocean 
structures. 
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I.B   BACKGROUND ; 

Attempts  to  analytically  predict  or  to  simulate  thermal 
stresses  and  metal  movement  which  occur  during  the  welding 
process  have  been  the  subject  of  considerable  research  effort 
during  the  past  ten  years.   The  complexity  of  the  interrela- 
tionships among  welding  parameters  which  contribute  to  the 
formation  of  residual  stresses  in  weldments  does  not  allow 
the  development  of  empirical  relationships  which  are  useful 
for  more  than  one  type  of  material,  joint  geometry,  or  welding 
process.   Welding  variables  were  too  numerous  and  their  in- 
terrelationships too  complex  for  the  development  of  a  general 
technique  for  the  analysis  of  transient  thermal  stresses, 
residual  stresses,  and  distortion  in  weldments w   Design 
engineers  usually  relied  upon  trial  and  error,  or  the  use  of 
design  lanes  derived  from  similar  applications  of  proposed 
welding  processes. 

Computer  simulation  of  the  welding  process  has  offered 
reasonable  opportunities  to  develop  techniques  which  apply 
to  wide  ranges  of  materials,  weld  geometries,  and  processes. 
A  description  of  the  chronological  and  systematic  development 
of  techniques  used  for  analysis  of  thermal  stresses  and  metal 
movement  during  the  welding  process  is  contained  in  a  paper 
presented  by  Masubuchi   to  an  international  conference  on 
Residual  Stresses  in  Welded  Construction  and  their  Effects. 
This  conference  was  held  in  London   15-17  November  1977, 
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hosted  by  the  Welding  Institute,  Abington  Hall,  Abington, 
Cambridge,  England.    Systematic  development  of  computer 

techniques  for  welding  simulation  began  with  a  one-dimen- 

2 

sional  analysis  by  Tall   in  1961.   This  research  involved 

the  simulation  of  thermal  stresses  occurring  during  the 

bead  welding  of  a  strip.   Efforts  at  Massachusetts  Institute 

3 
of  Technology  by  Masubuchi   in  1968  expanded  Tail's  analysis 

into  a  one-dimensional  FORTRAN  program.   This  research  was 

expanded  into  the  development  of  a  two-dimensional  FORTRAN 

4 
program  in  1970  by  Iwaki  .   Iwaki  attempted  to  apply  two- 
dimensional  finite  analysis  to  bead  on  plate  welding.   The 

most  recent  computer  program  to  be  used  at  MIT  was  developed 

5 
by  Muraki  .   Muraki ' s  program  uses  finite  element  analysis 

techniques  to  predict  thermal  stresses  and  metal  movement 

based  upon  temperature  distributions  within  the  weldment 

which  are  calculated  by  heat  flow  programs  also  developed 

at  MIT. 

I.C   SCOPE  AND  OBJECTIVES  OF  THIS  STUDY: 

The  purpose  of  this  research  is  to  design  and  carry 
out  an  experimental  procedure  to  generate  data  which  can 
be  used  to  validate  the  existing  heat  flow  program  and  the 
program  which  uses  finite  element  analysis  to  predict  the 
occurrence  of  thermal  stresses  and  distortion  during  welding 
The  production  of  seam  and  girth  welds  in  a  high  strength 
steel  cylinder  under  very  carefully  controlled  conditions 


17 

while  measuring  temperature,  strain,  and  distortion  as  a 
function  of  time  and  location  with  respect  to  arc  travel 
would  serve  as  a  valuable  aid  in  evaluating  the  accuracy 
and  potential  usefulness  of  the  heat  flow  program  and  of 
Muraki ' s  axisymmetric  finite  element  program  for  analysis 
of  thermal  stress  and  metal  movement  during  welding.   This 
research  is  part  of  a  continuing  effort  at  MIT  to  validate 
computer  analysis  of  welding.   A  summary  of  past  experimental 
studies  of  thermal  stresses  and  metal  movement  during  welding 
was  presented  by  Masubuchi   and  is  reproduced  and  updated 
here  as  Table  I-A. 
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CHAPTER  II 

MATERIAL  CONSIDERATIONS 
II. A   GENERAL: 

The  Office  of  Naval  Research  is  sponsoring  an  ongoing 
program  to  investigate  thermal  stresses  and  distortion  in 
high  strength  steel  (HY-130)  structures.   HY-130  steel  is  a 
high  quality,  low  carbon,  quenched  and  tempered  steel  with 
a  yield  strength  of  130-150  ksi.   This  material  was  produced 
by  the  United  States  Steel  Company,  primarily  for  use  by  the 
U.S.  Navy . 
II. B   COMPOSITION: 

Acceptable  ranges  of  chemical  composition  for  HY-13  0 

steel  are  shown  in  Table  II-A.   These  ranges,  reported  by 

20 
Shackleton,   are  compositional  limits  of  HY-130  steel  as 

received  for  production. 

Mechanical  properties  of  HY-130  steel  in  the  "as  received" 

21 

quenched  and  tempered  condition,  as  reported  by  Schrodt,   are: 

Yield  Strength:       130-145  ksi 

Elongation  in  2  inches:       15% 

Reduction  of  Area:        50%  minimum  transverse 

70%  through  thickness 

Impact  requirements:     60  ft-lbf  at  0   F 

A  Rockwell  hardness  test  conducted  on  the  material  to 

be  used  in  this  experimental  study  yielded  a  Rockwell  "C" 

Scale  Hardness  of  32. 
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COMPOSITIONAL  RANGES  OF  HY-130  STEEL  (%)20 


C  0.08  -  0.12 

Mn  0.60  -  0.90 

Si  0.20  -  0.35 

Ni  4.75  -  5.25 

Cr  0.40  -  0.70 

Mo  0.30  -  0.65 

V  0.05-0.10 

S  0.015  Max 

P  0.010  Max 

Ti  0.020  Max 

Cu  0.25  Max 


TABLE  I I -A 
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II. C   TEMPERATURE  DEPENDENCE  OF  MATERIAL  PROPERTIES; 

Subroutines  within  the  axi symmetric  finite  element 
program  developed  by  Muraki  for  analysis  of  thermal  stresses 
and  metal  movement  during  welding  require  that  material  pro- 
perties be  specified  as  a  function  of  temperature.  A  search 
of  the  literature  failed  to  produce  detailed  temperature 
relationships  among  the  various  material  properties  required 
by  the  Muraki  program  for  HY-130  at  elevated  temperatures. 

Graphical  presentation  of  these  relationships  reported  by 

21 
Schrodt   are  reproduced  in  Figures  II-l  through  11-10.   A 

tabular  summary  of  this  information  is  included  as  Tables  II-B 

and  II-C. 

The  Muraki  program  which  is  to  be  verified  by  experi- 
mental data  assumes  a  linear  relationship  between  various 
material  properties  and  temperature.   Assumed  linear  relation- 
ships which  served  as  input  data  to  the  Muraki  program  are 
indicated  by  a  dashed  line  on  the  appropriate  figures. 
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FIGURE  II-l   TYPICAL  STRESS-STRAIN  DIAGRAM  FOR 
HY-130  AT  ROOM  TEMPERATURE 
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FIGURE  II-2   ESTIMATED  EFFECT  OF  TEMPERATURE  ON  0.2' 
OFFSET  YIELD  STRESS  FOR  HY-130 
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FIGURE  II- 3   ESTIMATED  EFFECT  OF  TEMPERATURE  ON 
YOUNG'S  MODULUS  FOR  HY-130 
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FIGURE  II-4   TANGENT  MODULUS  CURVES  FOR  HY-130 
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FIGURE  II-5 


ESTIMATED  EFFECT  OF  TEMPERATURE  ON 
TANGENT  MODULUS  FOR  HY-130 
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FIGURE  II-6   ESTIMATED  EFFECT  OF  TEMPERATURE  ON 
POISSON'S  RATIO  FOR  HY-130 
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FIGURE  II-7 


ESTIMATED 
EXPANSION 


COEFFICIENT 
FOR  HY-130 


OF  THERMAL 


30 


en 


c 

•H 


>< 

Eh 

H 

O 

W 

o 

z 

CO 

w 

CN 

Q 

• 

o 
r- 

CM 


O 

o 
r- 


600 


1200  1800 

TEMPERATURE  (°F) 


2400 


FIGURE  II-8   EFFECT  OF  TEMPERATURE  ON  DENSITY  OF  HY-130 
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FIGURE  II-9 


ESTIMATED  EFFECT  OF  TEMPERATURE  ON 
THERMAL  CONDUCTIVITY  OF  HY-130 
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FIGURE  11-10 


ESTIMATED  EFFECT  OF  TEMPERATURE  ON 
SPECIFIC  HEAT  OF  HY-13  0 
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SUMMARY  OF  MECHANICAL  PROPERTIES 


FOR  HY-130  STEEL 


TEMPERATURE 
(°F) 

a 

E 

H 

V 

68 

143 

30 

22.2 

.30 

200 

138 

29.6 

21.8 

.308 

400 

131.5 

28.8 

21.2 

.321 

600 

122 

27.7 

20.5 

.334 

800 

110 

26 

19.2 

.346 

1000 

88 

22.7 

17.0 

.359 

1200 

50 

14.0 

9.2 

.372 

1400 

25,5 

5.0 

4.9 

.384 

1600 

11.2 

4.0 

3.1 

.397 

1800 

6.0 

3.0 

2.2 

.410 

2000 

4.0 

2.0 

1.5 

.422 

2200 

2.0 

1.2 

.8 

.435 

2400 

1.0 

.4 

.4 

.447 

2500 

0 

0 

0 

.454 

a  =   Yield  Stress  2%  Offset  (ksi) 

-3 
E  =   Elastic  Modulus  (ksi  x  10   ) 

-3 
H  =   Tangent  Modulus  (ksi  x  10   ) 

v  =   Poisson's  Ratio 


TABLE  II-B 
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SUMMARY  OF  PHYSICAL  PROPERTIES 


FOR  HY-130  STEEL 


TEMPERATURE 
(OF) 

P 

aT 

C 

K 

68 

0.284 

6.1 

19.4 

0.107 

200 

0.283 

6.5 

20.4 

0.11 

400 

0.282 

6.95 

21.5 

0.118 

600 

0.280 

7.3 

21.8 

0.126 

800 

0.279 

7.7 

21.2 

0.136 

1000 

0.278 

8.0 

19.5 

0.148 

1200 

0.276 

8.3 

17.1 

0.168 

1400 

0.275 

8.18 

15.0 

0.197 

1600 

0.275 

7.05 

13.4 

0.255 

1660 

0.275 

6.92 

13.1 

0.293 

1800 

0.275 

7.0 

14.0 

0.18 

2000 

0.275 

7.2 

15.0 

0.158 

2200 

0.275 

7.3 

16.0 

0.158 

2400 

0.275 

7.4 

16.7 

0.158 

2600 

0.275 

7.5 

17.0 

0.158 

2800 

0.275 

7.5 

7.0 

0.18 

p  =  Density  (lbm/in  ) 


a  =  Coefficient  of  Thermal  Expansion 

(micros train/o) 

C  =  Thermal  Conductivity 

(BTU/hr  ft  °F) 

k  =  Specific  Heat 

(BTU/lbm  °F) 

TABLE  II-C 
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CHAPTER  III 

DESCRIPTION  OF  EXPERIMENTAL  PROCEDURE 

III. A   FABRICATION  OF  THE  WELD  SPECIMEN; 

A  circular  cylinder  of  high  strength  steel  was  formed 
in  order  to  test  the  MIT  computer  program  for  the  analysis  of 
heat  flow  in  a  weldment,  with  subsequent  use  of  the  tempera- 
ture distributions  from  the  heat  flow  program  to  validate  the 
Muraki  axisymmetric  finite  element  program  for  analysis  of 
thermal  stress  and  metal  movement  during  welding.   The  HY-130 
circular  cylinder  was  formed  to  dimensions  suggested  in  the 
June  1977  proposal  on  Study  of  Residual  Stresses  and  Distor- 
tion in  Structural  Weldments  in  High  Strength  Steels.   This 
research  proposal  was  made  from  Massachusetts  Institute  of 
Technology  to  the  Office  of  Naval  Research. 

Material  for  forming  the  cylinder  was  obtained  from  the 
Naval  Ship  Research  and  Development  Center,  Annapolis  Labor- 
atory, Annapolis,  Maryland.   The  material  was  received  at 
Massachusetts  Institute  of  Technology  Welding  Laboratory  in 
the  form  of  HY-130  plates  60"  x  63"  x  3/4".   A  request  was 
made  to  Portsmouth  Naval  Shipyard  in  Portsmouth,  New  Hampshire, 
to  roll  the  HY-130  plates  into  cylindrical  shapes;  Portsmouth 
Naval  Shipyard  being  the  only  industrial  activity  in  the 
vicinity  of  MIT  with  experience  in  forming  HY-130.   Detailed 
drawings  and  specific  requirements  were  submitted  to  Portsmouth 
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Naval  Shipyard  and  an  agreement  was  reached  for  the  shipyard 
to  form  the  plates  into  cylinders.   However,  due  to  the 
thickness  (3/4")  of  the  plate  to  be  formed  to  a  small  outside 
diameter  (18"),  Portsmouth  Naval  Shipyard  engineers  suggested 
that  it -would  be  necessary  to  form  a  cylinder  by  bending  two 
plates  to  semi-circular  cylindrical  shapes  and  then  welding 
the  two  halves  together.   In  view  of  this  requirement  and  in 
order  to  facilitate  easier  handling  of  the  HY-130  plates,  the 
two  plates  were  flame  cut  to  four  plates  60"  x  3iy  x  3/4" 
at  the  MIT  Welding  Laboratory  prior  to  transporting  to 
Portsmouth,  New  Hampshire. 

Stages  of  cylinder  fabrication  are  depicted  in  Figures 
III-l  and  III-2.   In  order  to  prevent  heat  resulting  from 

the  flame  cutting  operation  from  affecting  thermal  properties 
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of  the  steel  close  to  the  proposed  weld,  Tsai    suggested  that 

flats  required  for  bending  be  removed  by  flame  cutting  no 
closer  than  V  to  the  proposed  leading  edges  of  the  semi- 
circular cylindrical  mating  surfaces.   The  flats  were  removed 
in  this  manner  and  remaining  material  to  be  removed  in  order 
to  form  a  60   Vee  groove  with  1/8"  root,  as  illustrated  in 
Figure  III-2,  was  removed  by  mechanical  grinding.   The  semi- 
circular cylinders  were  then  transported  back  to  the  MIT 
Welding  Laboratory  for  experimental  welding. 
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STAGES  OF  CYLINDER  FABRICATION 


Stage  1    -    HI-130  plate  with  dimensions  shown. 


Stage  2     -     Plate  flame  cut  to  dimensions  shown 

FIGURE    III-l 


Stage  h    - 


Stags  $     - 


Plate  halves  rolled  to  semi- 
circular cylindrical  shapes 
with  flats  approximately  1  3/li" 
remaining  on  leading  edges. 


Flats  removed  by  flame  cutting 
no  closer  than  1/U"  to  mating 
edges  of  semi-cylinders. 


lA"  excess  material  removed 
from  mating  edges  of  semi- 
cylinders  and  30°  angle  pro- 
duced by  mechanical  grinding. 


FIGURE    III-2 
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III.B   SELECTION  OF  WELDING  CONSUMABLES: 

The  type  of  welding  suggested  in  the  research  proposal 
for  girth  welding  of  cylindrical  shells  was  multi-pass  gas 
metal  arc  welding  CGMAW) .   This  type  of  welding  requires  the 
use  of  shielding  gas  and  filler  wire,  both  of  which  can 
significantly  affect  weld  quality.   The  most  common  type  of 
shielding  gas  used  in  GMAW  of  low  carbon  steel  is  Argon  with 
25%  C02 .   However,  in  the  case  of  HY-130,  experience  has 

shown  that  the  use  of  Argon  with  2  5%  C0~  leads  to  inadequate 

20 
fracture  toughness  properties  in  the  weld  deposit.   Shackleton 

recommends  rich  Argon  with  less  than  5%  Oxygen  or  10%  CO-  to 

avoid  deterioration  in  the  fracture  toughness  properties  of 

the  weld  deposit  caused  by  Oxygen  pickup  from  the  shielding 

gas.   Shielding  gas  selected  and  purchased  for  this  experi- 


men 


t  was  Argon  with  1%  0~ .   This  gas  was  reported  by 
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Radzirainski   as  the  gas  most  effective  in  shielding  HY-130 

welds. 

Filler  wire  selected  for  the  HY-130  cylindrical  seam 
and  girth  welds  was  1/16"  diameter  Linde  140S  wire.   This 
wire  was  developed  by  Union  Carbide  specifically  for  HY-130 
weldments  and  was  unanimously  endorsed  in  the  literature.   Dur- 
ing periods  when  the  filler  wire  was  not  in  use,  it  was  stored 
in  a  drying  oven  containing  desiccative  material.   This  step 
was  viewed  as  necessary  in  order  to  prevent  moisture  or  other 
impurities,  or  oxides,  from  collecting  on  the  wire  surface. 
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III.C   SELECTION  OF  WELDING  PARAMETERS: 

1.  The  variations  in  weldment  geometry  which  could 
be  considered  for  this  experiment  were  somewhat  restricted 
by  the  fact  that  cylinder  size  required  that  the  weld  be 

made  only  from  the  outside.   Based  upon  several  previous 
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welding  experiments  reported  by  Doty   and  the  above  men- 
tioned restriction,  a  single  60   Vee  groove  was  selected  as 
near  optimum  for  producing  the  seam  and  girth  welds  of  the 
cylinder.   Weld  joint  geometry  for  this  experiment  is 
depicted  below  in  Figure  III-3. 


x       o      / 
\^    60°  ^/ 

\ included  ' 

\  angle  / 


8V  inside  diameter 
1/8"  mating  surface 


FIGURE  III-3   WELD  JOINT  GEOMETRY  FOR  PRODUCING 
SEAM  WELDS  ON  HY-130  CYLINDER 
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2.  For  the  purposes  of  this  study,  heat  input  (Q)  is 

defined  in  the  conventional  way  as  the  product  of  arc  voltage 
and  welding  current  divided  by  arc  travel  speed,  as  expressed 
below  in  Equation  3-1: 

HEAT  INPUT  =  Q  =  V  X  I  (Equation  3-1) 

There  were  some  variations  in  welding  literature  regard- 
ing the  upper  level  of  heat  input  which  would  produce  a  good 
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quality  weld  in  3/4"  HY-130  steel.   Shackleton   recommended 

a  maximum  heat  input  of  40,000  Joules/inch.   This  upper  limit 
was  also  recommended  by  Nunez    and  Conner.    Radziminski 
recommended  a  maximum  of  45,000  Joules/inch.   Although  there 
was  some  variation  in  the  suggested  upper  limits  of  heat 
input,  there  seemed  to  be  general  agreement  that  4  0,000  Joules/ 
inch  was  a  reasonable  value  for  the  production  of  good  quality 
welds  and  this  value  was  selected  for  production  of  the  seam 
and  girth  welds  on  the  HY-130  cylinder.   Maximum  and  minimum 
allowable  values  of  pre-heat  and  interpass  temperature  were 
also  discussed  in  the  above  quoted  references;  general  recom- 
mendations were  more  varied,  but  based  on  a  rough  average  of 
recommended  values,  a  common  pre-heat  and  interpass  tempera- 
ture of  200  °F  was  selected, 

3.  Filler  wire  feed  rate  is  automatically  controlled 
by  the  GMAW  welding  machine.   A  wire  feed  rate  of  approxi- 
mately 87  inches/minute  of  the  1/16"  Linde  140S  filler  wire 
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was  measured  during  the  production  of  a  seam  weld.   The  flow 

3 
meter  for  the  Argon/1%  0      shielding  gas  was  set  at  40  ft  /hr. 

A  summary  of  welding  parameters  which  remained  unchanged  dur- 
ing the  course  of  this  experiment  is  shown  in  Table  III-A 
below.   The  only  parameters  which  were  varied  between  welding 
passes  were  the  arc  voltage,  welding  current,  and  in  one 
instance  the  arc  travel  speed.   Voltage,  current,  and  arc 
travel  speed  were  recorded  for  each  pass  and  are  included  in 
the  presentation  of  results. 


WELDING  PARAMETERS 


JOINT  GEOMETRY  60°  Single  V 

Groove 

MAXIMUM  HEAT  INPUT  40,000  Joules/in 

FILLER  WIRE  1/16 "  Diameter 

Linde  140S 

feed  rate  -  87  in/min 

SHIELDING  GAS  Argon/1%  02 

flow  rate  -  40  ft3/hr 


TABLE  III-A 


III.D   ARRANGEMENT  OF  APPARATUS: 

Welding  equipment  used  for  the  entire  experiment 
consisted  of  a  Linde  SVI-300  Power  Supply  with  ST-2  Torch 
and  associated  arc  travel  and  wire  feed  mechanisms.   Seam 
welds  were  produced  by  placing  the  cylinder  on  an  adjustable 
welding  table  and  moving  the  welding  torch  along  the  groove 
of  the  seam  weld.   Prior  to  positioning  the  cylinder  on  the 
welding  table,  the  two  cylinder  halves  were  tack  welded  at 
positions  6"  from  either  end  and  at  the  center  with  1"  welds 
on  the  inside  of  each  mating  edge.   Positions  of  the  cylinder 
and  welding  torch  for  production  of  a  seam  weld  are  depicted 
in  the  photograph,  Figure  III-4,  below. 


FIGURE  III-4   POSITION  OF  CYLINDER  AND  WELDING  TORCH  FOR 
PRODUCTION  OF  SEAM  WELD 
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III.E   TESTING  OF  WELDING  EQUIPMENT; 

In  order  to  determine  approximate  temperature  distri- 
butions during  actual  welding  and  to  facilitate  a  thorough 
checkout  of  welding  equipment  prior  to  welding  passes  for 
data  gathering,  it  was  deemed  advantageous  to  produce  a  seam 
weld  on  one  side  of  the  cylinder  without  benefit  of  instru- 
mentation to  record  temperature  and  strain.   This  allowed 
the  wire  feed  mechanism,  arc  travel  controls  and  indicators, 
gas  flow  regulators  and  meters,  and  the  welding  equipment 
electrical  circuity  and  instrumentation  to  be  tested  and 
checked  for  accuracy  under  conditions  expected  during  the 
actual  test.   This  "test  run"  also  provided  an  opportunity 
for  operator  training.   A  preliminary  checkout  of  welding 
equipment  had  been  conducted  using  "bead-on-plate"  welds  on 
scrap  pieces  of  HY-130. 

Approximate  temperature  distributions  to  be  expected 
during  the  instrumented  test  were  determined  by  the  applica- 
tion of  temperature  indicating  liquid  (TEMPILAQ)  in  narrow 
bands  transverse  to  the  welding  direction.   TEMPILAQ  liquids 
corresponding  to  temperature  levels  of  400   F,  600   F,  800   F, 
1000   F,  and  1100   F  were  applied  to  the  chemically  cleaned 
inner  and  outer  surfaces  of  the  cylinder.   TEMPILAQ  was  re- 
applied between  passes  and  distances  of  temperature  penetra- 
tion from  the  weld  centerline  on  both  the  inside  and  outside 
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of  the  cylinder  were  recorded  for  each  welding  pass.   Welding 
current,  arc  voltage,  arc  travel  speed,  and  maximum  distances 
from  the  weld  centerline  for  the  various  TEMPILAQ  levels  were 
recorded  for  each  welding  pass  and  are  shown  in  Table  III-B. 
This  tabulated  data  was  used  as  a  basis  for  determining  strain 
gage  and  thermocouple  locations  for  the  instrumented  seam  and 
girth  welds  to  be  used  in  validating  heat  flow  and  thermal 
stress  computer  programs. 
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WEEDING  PASS 

♦  TEMPERATURE 

PENETRATION  DISTANCE 

OUTSIDE 

INSIDE 

LEFT   RIGHT 

LEFT 

RIGHT 

1.  I  =  260A 

UOOF      .86 

.9U 

.88 

7=  23V 

600F 

.63 

.60 

t  =  12  in/min 

800F 

.uu 

.UO 

1000F 

.38 

.35 

HOOF 

.25 

.22 

2.  I  =  260A 

UOOF     .61     .38 

.69 

.75 

V  ■  23V 

600F 

.38 

.56 

v  »  12  in/min 

800F 

.22 

.33 

1000F 

.06 

.31 

1100F 

- 

.19 

3.  I  =  260A 

UOOF     .80 

.88 

.63 

V  -  23V 

600F 

.50 

.38 

v  ■  12  in/min 

800F 

.23 

- 

1000F 

- 

- 

HOOF 

- 

- 

U.  I  a  26QA 

UOOF      -      .92 

.50 

.83 

V  =  23V 

600F      -      .U9 

.19 

.50 

t  =  12  in/min 

800F 

- 

- 

1000F 

- 

- 

1100F 

- 

- 

5.  I  =  350A 

UOOF     .92    1.55 

1.06 

1.63 

V  =  22V 

600F      -      .80 

.38 

.75 

v  =  10  in/min 

800F      -      .58 

- 

.38 

1000F      -      .09 

- 

- 

1100F      -      .03 

- 

- 

6.  I  -  330A 

UOOF     1.U9    1.15 

1.50 

.88 

V  =  21V 

600F     .30     .75 

.75 

- 

v  -  10  in/  min 

800F     .55 

.25 

— 

1000F     ,k9 

- 

- 

HOOF      .h9 

- 

- 

7.  I  =  350A 

UOOF     .U6     .78 

.38 

.75 

V  =  22V 

600F      -      .UO 

- 

.13 

v  »  20  in/min 

800F      -      .25 

- 

- 

1000F      -      .03 

- 

- 

HOOF 

- 

- 

8.  I  =  370A 

TEMPERATURE  PENETRATION  DISTANCES  NOT 

RECORDED 

V  =  21V 

v  »  20  in/min 

*  Temperature  pe  net  rat 

,ion  distances  are  expressed 

in  inches, 

measured 

in  the  direction  inc 

Licated  from  the  weld  centerline. 

TABLE    III-B 
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III.F   SELECTION  AND  ARRANGEMENT  OF  INSTRUMENTATION  TO 
MEASURE  STRESS  AND  DISTORTION  DURING  WELDING: 

Stresses  occurring  during  the  thermal  cycles  associated 
with  seam  and  girth  welding  of  the  HY-130  cylinder  were  de- 
termined by  measuring  strains  as  a  function  of  time  (i.e.  - 
location  of  the  strain  measurement  with  respect  to  the  time 
history  of  welding  arc  location)  and  as  a  function  of  distance 
from  the  weld  centerline.   Strains  were  measured  by  placing 
electrical-resistance  wire  gages  on  the  outside  surface  of 
the  cylinder.   Bi-axial  strain  gages  were  placed  with  one 
axis  aligned  with  the  cylinder  axis  and  welding  direction 
(herein  called  the  longitudinal  direction)  and  the  other  axis 
aligned  in  a  direction  normal  to  the  welding  direction  (herein 
referred  to  as  the  transverse  direction.   Characteristics  of 
the  strain  gages  used  for  this  experiment  are  shown  in 
Table  III-C. 

In  order  to  validate  the  MIT  heat  flow  program,  it  was 
necessary  to  measure  the  heat  distribution  within  the  cylinder 
Heat  flow  and  temperature  distribution  within  the  HY-130 
cylinder  was  viewed  as  critical  to  validation  of  the  Muraki 
program,  inasmuch  as  strain,  stress,  and  distortion  calcula- 
tions are  all  based  upon  temperature  distribution  data  which 
must  be  supplied  as  input  to  the  program.   Thermocouples  used 
for  measuring  temperature  distribution  were  Chromel/Alumel 
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STRAIN  GAGE  CHARACTERISTICS 

Type  Used: 

FAET-12D-12S6(EL) 

Manufacturer : 

BLH  Electronics,  Inc. 
Waltham,  Massachusetts  02154 

Gage  Factor: 

2.00  i  .1% 

Resistance : 

120  +  .2  Ohms 

Grid  Length: 

0.125  inches 

Cement  Used: 

BLH  EPY-500 

TABLE  III-C 

Strain  indications  were  corrected  for  temperature 
by  the  following  temperature  error  (.apparent  strain)  equation: 


EAP  =  -91.37  +  2.17T  -  2.33  x  10~2T2  +  5.31  x  10  5T3 

-8  4 
-  2.40  x  10   T 
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wires  which  were  peened  into  small  (0.040")  indentations 
drilled  into  the  cylinder  outer  surface.   Thermocouple 
circuit  continuity  was  assured  by  the  use  of  an  inductive 
spot  welder  on  the  thermocouple  connections.   After  circuit 
continuity  was  assured,  the  thermocouples  were  protected  by 
the  use  of  sealing  cement. 

Since  longitudinal  and  transverse  strains  were  to  be 
measured  with  gages  which  required  temperature  correction, 
it  was  necessary  to  place  thermocouples  in  positions  which 
would  give  temperature-time  histories  consistent  with  strain 
gage  locations.   Arrangement  of  strain  gages  and  thermo- 
couples on  the  HY-130  cylinder  outer  surface  is  shown  in 
Figure  III-5.   Arranging  the  strain  gages  and  thermocouples 
in  this  manner  allowed  them  to  be  used  for  measuring  strain 
and  temperature  distributions  during  the  production  of  the 
second  seam  weld  (the  first  seam  weld  was  made  without  in- 
strumentation) and  then  to  be  used  again  after  the  cylinder 
had  been  cut  into  sections  and  mating  edges  prepared  for  a 
girth  weld.   Arranging  the  strain  gages  and  thermocouples 
at  a  45   angle  with  respect  to  the  seam  weld  and  at  equi- 
distant locations  with  respect  to  each  other  allowed  the 
strain  gages  and  thermocouples  to  give  strain-temperature 
distributions  which  could  be  adjusted  in  time  so  as  to  be 
coincident. 
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ARRANGEMENT  OF  STRAIN  GAGES  AND  THERMOCOUPLES 


Z^\ 


seam 


u  c  c  mm  m  I  <  ( ( c  g  ES33& 


FIGURE    III-5 
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It  should  be  pointed  out  that  the  practice  of  locating 
strain  gages  and  thermocouples  so  as  to  require  time  adjust- 
ment of  their  readings  assumes  quasi-steady-state  heat  flow 
(i.e.  -  the  heat  flow/temperature  distribution  with  respect 
to  the  welding  arc  is  the  same  at  a  position  40"  along  the 
length  of  the  cylinder  as  at  a  position  41"  along  the  cylin- 
der) .   This  assumption  is  necessary  in  view  of  the  fact  that 
strain  gages  and  thermocouples  cannot  be  located  at  the  same 
position. 

Arranging  the  strain  gages  and  thermocouples  at  a 
45   angle  with  respect  to  the  seam  weld  was  viewed  as  ad- 
vantageous since  this  scheme  places  the  strain  gages  and  t 
thermocouples  which  were  subjected  to  the  most  severe  thermal 
cycle  during  production  of  the  seam  weld  in  a  position  which 
is  subjected  to  the  least  severe  thermal  cycle  during  the 
girth  weld.   Failures  occurring  during  the  seam  weld  are 
most  likely  to  occur  at  locations  which  are  least  critical 
to  measurements  to  be  taken  during  the  girth  weld. 

After  determining  and  marking  the  desired  strain  gage 
locations,  strain  gages  were  cemented  to  the  cylinder  outer 
surface.   Physical  size  and  weight  of  the  cylinder  prevented 
the  usual  technique  of  curing  the  strain  gage  bonding  cement 
in  laboratory  ovens.   In  order  to  assure  that  the  strain  gage 
bonding  would  withstand  expected  temperatures,  a  5  foot  by 
5  foot  thermal  insulating  blanket  (graciously  manufactured 
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for  this  purpose  by  Bethlehem  Steel  Corporation  Shipyard  in 
Boston,  Massachusetts)  was  placed  around  the  cylinder  and 
electrical  heating  elements  placed  inside.   The  temperature 
of  the  cylinder  was  slowly  raised  to  250   F  and  the  strain 
gage  bonding  cement  was  cured  at  this  temperature  for  a  period 
of  four  hours,  as  directed  by  the  strain  gage  manufacturer's 
instructions.   The  heating  elements  used  for  this  purpose 
were  later  used  to  preheat  the  cylinder  prior  to  welding. 

For  production  of  the  girth  weld,  it  was  necessary  to 
cut  the  HY-130  cylinder  into  sections.   Original  plans  called 
for  cutting  the  60  inch  cylinder  into  two  30  inch  sections 
which  would  have  mating  edges  prepared  for  a  girth  weld  by 
machining  a  60   single  Vee  groove.   Due  to  physical  limitations 
of  laboratory  facilities  it  was  necessary  to  reduce  the  length 
of  the  two  sections  to  be  joined  by  a  girth  weld  to  20  inches 
each.   This  required  that  the  HY-13  0  cylinder  be  cut  into 
three  sections,  each  20  inches  in  length.   This  decision  was 
made  prior  to  production  of  the  instrumented  seam  weld,  and 
the  strain  gages  were  placed  so  as  to  permit  measurements 
during  the  seam  weld  as  well  as  during  production  of  the  girth 
weld. 

After  completion  of  the  seam  weld,  the  60  inch  long 
HY-130  cylinder  was  transported  to  the  Army  Materials  and 
Mechanics  Research  Center  in  Watertown,  Massachusetts,  where 
it  was  sawed  into  three  sections,  each  2  0  inches  in  length. 
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A  large  power  hack  saw  was  used  to  cut  the  cylinder  into 
sections.   The  saw  was  used  on  slow  speed  and  the  blade  and 
cylinder  were  continuously  cooled  by  an  oil  and  water  cool- 
ing mixture  during  cutting.   This  procedure  was  viewed  as 
necessary  in  order  to  prevent  heat  generated  by  the  cutting 
action  from  affecting  thermal  properties  of  the  mating  edges 
to  be  joined  by  girth  welding. 

After  the  two  saw  cuts  were  made,  producing  the  desired 
three  20  inch  long  cylindrical  sections,  a  %  inch  thick  ring 
was  sawed  from  the  section  which  would  not  be  used  in  the 
girth  weld;  again  paying  particular  attention  to  preventing 
heat  buildup  in  the  material  being  sawed.   A  diagram  of  the 
resulting  sections  appears  below  as  Figure  III-6. 


*  Approximately  h   inch  of  cylinder  length 
lost  to  width  of  saw  cuts 


FIGURE  III-6 
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A  cross-section  of  the  previously  produced  seam  weld 
was  removed  from  the  %   inch  thick  ring  and  later  ground, 
polished,  and  etched  for  analysis  of  the  metallurgical  proper- 
ties of  the  seam  weld.   Details  of  this  procedure  and  a  dis- 
cussion of  the  results  are  contained  in  Chapter  IV  of  this 
study. 

Mating  edges  of  the  two  sections  to  be  joined  by  girth 
welding  were  then  prepared  for  welding  by  machining  a  60° 
single  Vee  groove  as  shown  below  in  Figure  III-7.   Heat 
buildup  during  machining  was  prevented  by  the  use  of  a  con- 
tinuously flowing  oil  and  water  cooling  bath  during  machining, 


20  in.- 


to  cylinder  end 


1/8  inch 


2  0  in. 


to  cylinder  end 


FIGURE  III-7 
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The  two  20  inch  cylindrical  sections  were  transported 

back  to  the  MIT  welding  laboratory  where  apparatus  required 
for  production  of  the  girth  weld  was  assembled.   In  order  to 
control  arc  travel  speed  during  production  of  the  girth  weld, 
it  was  necessary  to  design  and  manufacture  a  supporting  struc- 
ture for  the  cylinder  which  would  permit  cylinder  rotation  at 
a  controlled  speed  with  respect  to  the  fixed  position  of  the 
GMAW  torch.   This  was  accomplished  by  producing  a  structure 
which  was  designed  to  support  the  two  20  inch  cylinders  from 
within.   Details  of  this  supporting  structure  are  shown  below 
in  Figure  III-8. 


FIGURE  III-8 
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The  rectangular  base  of  the  supporting  structure  was 
bolted  to  the  base  of  a  tilting  turntable  and  the  20  inch 
long  cylinder  halves  were  fixed  to  the  support  by  force 
exerted  on  the  interior  walls  of  the  cylinder  by  the  threaded 
dowels.   The  turntable  was  then  tilted  so  as  to  place  the 
now  40  inch  long  cylinder  (with  60°  Vee  groove  located  in  the 
center)  in  a  horizontal  position.   The  threaded  dowels  were 
adjusted  in  order  to  provide  for  constant  clearance  between 
the  stationary  GMAW  torch  and  cylinder  outer  wall  during 
cylinder  rotation.   Turntable  controls  were  adjusted  and 
calibrated  in  order  to  achieve  desired  arc  travel  speeds. 
A  photograph  of  the  apparatus  in  use  during  production  of  the 
girth  weld  appears  below  as  Figure  III-9. 


FIGURE  III-9 
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In  order  to  measure  radial  distortion  during  production 

of  the  girth  weld,  two  Bourns  Linear  Motion  Potentiometers 
were  placed  at  the  inner  surface  of  the  cylinder.   Character- 
istics of  the  linear  motion  potentiometers  are  shown  below 
in  Table  III-D. 


ACCEPTANCE  TEST 

RECORD 

FOR  LINEAR  MOTION 

POTENTIOMETERS 

TESTS  PERFORMED 

SPEC  LIMITS 

MECHANICAL  TRAVEL  

0 . 44  in  max 

ACTUATION  FORCE  

4  oz  min 

TOTAL  ELEMENT  RESISTANCE 
INSULATION  RESISTANCE 

.  .  4750  to  5250  Ohms 

50  Megohms  min 

at  500  V  DC 

NOISE  OR  CONTACT 

RESISTANCE  VARIATION.. 

10  0  ohms  max 

END  SETTINGS  

0.1  to  5.0% 
each  end 

LINEARITY  

*  1.0%  VR  max 

TABLE  III-D 
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Displacement  transducer  circuitry  and  arrangement  of 
linear  motion  potentiometers  within  the  HY-130  cylinder  are 
indicated  in  Figures  111-10  and  III-ll  respectively.   Prior 
to  mounting  the  displacement  transducers  in  position  within 
the  cylinder,  a  calibration  test  was  conducted  by  deflecting 
the  displacement  transducers  through  measured  distances  and 
recording  resultant  voltages.   A  plot  of  this  calibration 
test  appears  as  Figure  111-12.   After  completion  of  calibra- 
tion tests,  the  linear  motion  potentiometers  were  installed 
within  the  cylinder  as  indicated  in  Figure  III-ll.   Holes 
were  drilled  and  threaded  in  the  cylinder  support  mechanism 
and  the  linear  motion  potentiometer  support  rods  were 
threaded  and  screwed  into  these  holes  at  diametrically 
opposed  locations  within  the  cylinder.   Potentiometers 
corresponding  to  Channels  1  and  2  of  the  recorder  were 
located  at  positions  1/2  inch  and  1  inch  respectively  in 
the  cylinder  axial  direction  from  the  weld  centerline. 
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6 


6 


6- 


**- 


6V  LU 

DC  ' 


It 


GOULD  BRUSH  220 
TWO-CHANNEL  RECORDER 


FIGURE  111-10   DISPLACEMENT  TRANSDUCER  CIRCUITRY 


CHANNEL  1  DISPLACEMENT  TRANSDUCER 
SENSING  ELEMENT  LOCATED  H    INCH 
FROM  WELD  CENTERLINE.   INSULATED 

FROM  CYLINDER  INNER  SURFACE 

BY  CERAMIC  BEAD. 


DISPLACEMENT  TRANSDUCER  SENSING  E 
MOUNTING  RODS,  SCREWED  INTO 
CYLINDER  SUPPORT. 


CHANNEL  2  DISPLACEMENT  TRANSDUCER 
SENSING  ELEMENT  LOCATED  1  INCH  FROM 
WELD  CENTERLINE.   ALSO  INSULATED. 


CYLINDER  SUPPORT  MECHANISM 


FIGURE  III-ll 


DISPLACEMENT  TRANSDUCER  SENSING  ELEMENT 
LOCATIONS 
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Strain  and  temperature  were  measured  at  the  indicated 

locations  during  production  of  the  instrumented  seam  weld; 

and  strain,  temperature,  and  distortion  were  measured  during 

production  of  the  girth  weld.   Equipment  used  to  record  the 

measured  values  included  the  following: 

Minneapolis-Honeywell  H1508  24  Channel  Visicorder 

Three  Helland  Universal  Bridge  Balancer  Units, 
Model  no.  82-6  (ser  no's  1332,1307,1312) 

16  Channel  Strain  Gage  Calibrator 
(one  arm  bridge  adapter) 

Gould  Brush  220  Two-Channel  Recorder  (distortion) 

General  arrangement  and  operation  of  this  equipment 
during  production  of  the  girth  weld  are  indicated  below  in 
Figure  111-13. 


FIGURE  111-13 
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III.G   PRODUCTION  OF  THE  INSTRUMENTED  SEAM  AND  GIRTH  WELDS: 

1.   The  Seam  Weld:   After  the  HY-130  cylinder  had  been 
suitably  supported  and  aligned;  instrumentation  and  other 
apparatus  arranged,  checked  and  calibrated,  the  instrumented 
seam  weld  was  produced  by  making  seven  passes.   Figure  111-14 
below  shows  approximate  dimensions  and  location  of  each  pass. 
Values  of  weld  parameters  for  each  pass  are  tabulated  in 
Table  III-E. 


FIGURE  111-14 
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2.   The  Girth  Weld:   Production  of  the  girth  weld 
required  careful  positioning  of  the  rotating  HY-130  cylinder 
with  respect  to  the  GMA  welding  torch.   The  girth  weld  required 
six  passes  with  approximate  dimensions  shown  below  in  Figure 
111-15.   Girth  weld  parameters  are  also  included  in  Table 
III-E. 


FIGURE  111-15 
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SEAM  WELD  PASSES 

GIRTH  WELD  PASSES 

WELD 
PASS  NO. 

I 

V 

V 

I 

V 

V 

1. 

320 

26 

0.20 

270 

26 

0.20 

2. 

320 

25 

0.20 

290 

25.5 

0.20 

3. 

320 

25 

0.20 

270 

26 

0.20 

4. 

300 

26 

0.27 

285 

25.5 

0.20 

5. 

285 

25 

0.20 

280 

26 

0.20 

6. 

280 

25 

0.20 

290 

26 

0.20 

7. 

290 

25 

0.40 

- 

- 

- 

I  =  Welding  Current  in  Amperes 

V  =  arc  voltage  in  Volts 

v  =  arc  travel  speed  in  inches/sec 


TABLE  III-E   SEAM  AND  GIRTH  WELD  PASSES 
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CHAPTER  IV 

METALLURGICAL  CONSIDERATIONS 
IV. A   PREPARATION  OF  SPECIMEN: 

When  the  60  inch  cylinder  was  sawed  into  sections,  as 
previously  discussed  in  Chapter  III,  a  cross-section  of  the 
instrumented  weld  was  cut  from  the  h   inch  ring.   The  weld 
specimen  was  cut  to  approximately  3/4  inch  by  3/4  inch  and 
marked  for  identification  as  shown  in  Figure  IV-1.   Mounting 
for  polishing  was  in  thermo-setting  bakelite  resin,  pressed 
into  a  V%   inch  diameter  by  3/4  inch  long  cylindrical  shape 
under  hydraulic  pressure  of  4200  psi  at  250   F  for  2  minutes 
using  a  Buehler  press. 

The  specimen  was  then  rough  polished  by  successive  use 
of  60,  120,  180,  240,  320,  and  400  grit  3M  TRI-M-ITE  Silicon 
Carbide  paper.   Fine  polishing  was  accomplished  by  successive 
use  of  l.Ou,  0.3y,  and  0.05y  alumina  polishing  powder.   Etching 
for  photo-micrographs  required  about  15  seconds  with  5%  NITAL. 
The  500X  photo-micrographs,  Figures  IV-2  through  IV- 5 ,  were 
produced  by  the  use  of  a  LEITZ  WETZLAR  microscope  with  NPL 
50X/0.85  lens,  using  Polaroid  type  52  POLAPAN  film.   Fine 
polishing  was  repeated  and  the  specimen  was  etched  using  1% 
NITAL  prior  to  the  production  of  1200X  and  6000X  electron  scan 
micrographs.   These  electron  scan  micrographs  appear  as  Figures 
IV-6  through  IV-13.   Locations  of  the  photomicrographs  and  the 
electron  scan  micrographs  are  indicated  on  Figure  IV-1. 
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©  - 
©  - 

©- 
©- 


Weld  Metal 


Fusion  Zone 


Heat  Affected  Zone 


Parent  Metal 


FIGURE  IV- 1 


DIAGRAM  SHOWING  LOCATIONS  OF  PHOTO-MICROGRAPHS 
AND  ELECTRON  SCAN  MICROGRAPHS.   3/4"  SAMPLE 
WAS  CUT  FROM  THE  SEAM  WELD  (SHOWING  CROSS 
SECTION  OF  WELD  -  WITH  NOTCH  IN  UPPER  LEFT 
CORRESPONDING  TO  WELD  CENTERLINE) .   AFTER 
POLISHING  AND  ETCHING, 50 OX  PHOTO-MICROGRAPHS 
AND  120 OX/6 00 OX  ELECTRON  SCAN  MICROGRAPHS  WERE 
PRODUCED  AT  THE  LOCATIONS  INDICATED. 
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FIG  IV- 2  PHOTOMICROGRAPH  OF  HY-130 

WELD  METAL.  50 OX.  5%  NITAL. 


FIG  IV- 3  PHOTOMICROGRAPH  OF  HY-130 

FUSION  ZONE.  500X.  5%  NITAL. 
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FIG  IV- 4  PHOTOMICROGRAPH  OF  HY-130 
HAZ.  500X.   5%  NITAL. 


■  ■     ■         '.  ---     '  \T.  .'' 

■  •  ;  ■  i       ^  . a. 


- 


* 
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FIG  IV- 5  PHOTOMICROGRAPH  OF  HY-130 

PARENT  METAL.  50 OX.  5%  NITAL. 
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FIG  IV- 6  ELECTRON  SCAN  MICROGRAPH  OF 
HY-130  WELD  METAL.   1200X. 
1%  NITAL. 


FIG  IV- 7  ELECTRON  SCAN  MICROGRAPH  OF 
HY-130  WELD  METAL.   60  00X. 
1%  NITAL. 
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FIG  IV- 8  ELECTRON  SCAN  MICROGRAPH  OF 
HY-130  FUSION  ZONE.  1200X. 
1%  NITAL. 


FIG  IV- 9  ELECTRON  SCAN  MICROGRAPH  OF 
HY-130  FUSION  ZONE.  6000X. 
1%  NITAL. 


FIG  IV-10  ELECTRON  SCAN  MICROGRAPH  OF 

HY-130  HAZ.   1200X.   1%  NITAL. 


FIG  IV-11  ELECTRON  SCAN  MICROGRAPH  OF 

HY-130  HAZ.   6000X.  1%  NITAL. 
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FIG  IV-12  ELECTRON  SCAN  MICROGRAPH  OF 
HY-130  PARENT  METAL.   1200X. 
1%  NITAL. 


FIG  IV- 13  ELECTRON  SCAN  MICROGRAPH  OF 
HY-130  PARENT  METAL.  6000X 
1%  NITAL. 
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CHAPTER  V 
DESCRIPTION  OF  ANALYTICAL  TOOLS  AND  PROCEDURES 

V.A   INTRODUCTION; 

Analytical  procedures  used  in  this  study  to  predict 
temperature  distribution,  strain,  stress,  and  distortion 
occurring  during  production  of  a  weldment  were  developed  at 
Massachusetts  Institute  of  Technology.   Both  the  heat  flow 
program  and  the  program  used  to  analyze  thermal  stress  and 
metal  movement  use  finite  element  techniques  and  are  written 
in  the  FORTRAN-IV  language. 

Discussion  of  these  programs  within  this  study  will  be 
limited  to  a  brief  synopsis  of  their  theoretical  background, 
presentation  of  assumptions  and  applicable  equations  upon 
which  they  are  based,  and  a  summary  of  welding  parameters 
pertinent  to  production  of  the  HY-130  weldment  which  were 
used  as  input  to  the  programs.   Detailed  discussion  of  these 
programs,  lists,  card  decks,  and  manuals  describing  their  use 
are  located  within  the  Department  of  Ocean  Engineering  at 
Massachusetts  Institute  of  Technology. 

V.B   DESCRIPTION  OF  THE  HEAT  FLOW  PROGRAM; 

1.   Theoretical  Background.   The  finite  element  computer 
program  used  in  this  study  to  generate  temperature  distribution 

in  the  HY-13  0  cylinder  girth  weld  is  Program  B-4  in  a  report 

27 
developed  by  Muraki   and  Masubuchi  entitled  "Computer  Programs 
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Useful  for  the  Analysis  of  Heat  Flow  in  Weldments"  at  MIT  in 

1974.   This  program  develops  temperature  distribution  based 

on  heat  conduction  in  quasi-stationary  state  from  a  moving 

point  heat  source  in  a  plate  of  finite  thickness.   The 

solution  assumes: 

a.  The  size  of  the  heat  source  is  small  compared 
to  that  of  the  body  under  consideration. 

b.  Heat  is  supplied  instantaneously. 

c.  Thermophysical  properties  do  not  depend  upon 
temperature . 

d.  Heat  is  transferred  in  conduction  only  (other 
modes  such  as  radiation  and  convection  are  not  considered) . 

e.  The  material  of  the  body  is  uniform  and  iso- 
tropic in  the  thermophysical  properties. 

The  solution  is  based  upon  the  following  equation  which 
expresses  the  difference  between  the  temperature  at  time  t 
tin  seconds  after  the  arc  is  initiated)  and  the  initial 
temperature: 


-vg 

T-T  =  tt3-  e  2K 
o   2ttc 


-vr     -vr, 

e  2Kn  +  e  2K 

r      r ' 
n      n 


(Eq.  5-1) 


Where :  £  =  x  -  vt 

2  2     2     2 

r   =  x  +  y   +  z 

r   =  z2  +  y2  +  (2nh  -  z)2 

n  2     2             2 

r'  =  C  +  y  +  (2nh  +  z) 

q  =  heat  intensity 
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2 .   Geometry  of  Weld  Specimen  and  Description  of 
Finite  Element  Representation  of  Temperature  Distribution 
Program: 

Division  of  the  HY-13  0  cylinder  girth  weld  into  elements 
for  use  in  the  heat  conduction  program  is  illustrated  in 
Figure  V-l.   Node  coordinates  (x,  y,  and  z  values)  for  each 
of  the  114  nodes  are  listed  in  Table  V-A.   The  weldment  is 
represented  by  finite  elements  which  are  arranged  so  as  to 
cause  node  locations  for  analytical  determination  of  heat 
distribution  to  be  coincident  with  locations  of  thermocouples 
which  were  actually  placed  on  the  cylinder  outer  surface  in 
order  to  experimentally  determine  temperature  distribution. 

Node  coordinates  listed  in  Table  V-A  are  based  upon  an 
origin  Cx=y=z=0)  which  is  coincident  with  the  starting  point 
of  the  heat  source  (welding  arc) ,  with  subsequent  times  of 
interest  corresponding  to  location  of  the  heat  source  (in 
nodal  coordinates)  at  the  referenced  time.   This  analysis 
assumes  that  the  heat  source  moves  on  the  plane  z=0  along 
the  x-axis.   The  x  direction  (direction  of  welding  arc 
travel)  would  appear  on  Figure  V-l  as  normal  to  the  y-z 
plane,  positive  direction  into  the  paper.   In  this  figure 
the  y-axis  is  along  the  cylinder  length  and  the  z-axis  is 
in  the  radial  direction. 
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In  order  to  determine  temperature  distribution  at 
weldment  nodes  for  the  first  pass  of  the  HY-130  cylinder 
girth  weld,  the  following  parameters  were  used  as  input  to 
the  heat  flow  program: 

Velocity  of  welding  arc  (v)  =  0.2000  in/sec 

Heat  intensity  (q)  =  6084  Joules/ 

J      ^  '  sec  in 

Where:   q  =   T  X  V  X  narc  (Eq.  5-2) 

h 

I  =  welding  current 

V  =  arc  voltage 

h  =  plate  thickness 

n    =  welding  arc  efficiency 
arc 

Density  (p)  =  0.2  82  lbm/.  3 

Specific  heat  (k)  =  124.4664  Joules/,,   oF 

Thermal  conductivity  (c)  =  0.5250  Joule/       o_ 

2  '  sec  in   F 

It  should  be  noted  that  this  stage  of  the  computer 

analysis  does  not  provide  for  temperature  dependence  of 

material  properties .   Values  used  in  this  case  to  determine 

temperature  distribution  are  those  valid  at  400   F.   It  is 

also  significant  to  note  that  for  this  analysis,  an  arc 

efficiency  (n    )  of  0.65  was  assumed, 
arc 
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No  modifications  were  required  to  the  analytical  or 
mathematical  techniques  required  for  generating  temperature 
distributions  to  be  used  by  the  Muraki  strain,  stress,  and 
distortion  calculating  program.   However,  since  large  amounts 
of  numerical  data  were  generated  by  the  heat  flow  program  for 
input  to  the  Muraki  program,  it  was  necessary  to  modify  the 
heat  program  output.   Since  weldment  geometry  and  finite 
element  representation  of  the  nodes  and  elements  were  the 
same  (except  for  coordinates)  for  both  programs,  the  heat  flow 
program  was  modified  so  as  to  provide  its  output  in  the  form 
of  punched  cards  suitable  for  use  as  input  to  the  Muraki 
program.   This  reduced  the  probability  of  introducing  numeri- 
cal error  betv/een  the  two  programs  as  well  as  improving  the 
ease  of  program  interaction. 
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V.C   DESCRIPTION  OF  THE  PROGRAM  TO  CALCULATE  STRAIN,  STRESS, 
AND  DISTORTION: 

1.   Theoretical  Background:   The  program  used  in  this 
study  to  analytically  generate  strain,  stress,  and  distortion 
distributions  within  the  HY-130  cylinder  girth  weld  was  de- 
veloped at  MIT  by  Toyohiki  Muraki  for  thermo  elasto-plastic 
analyses  of  axisymmetric  structural  components.   The  Muraki 
program  is  written  in  FORTRAN  IV  language  and  uses  finite 
element  techniques  in  an  incremental  fashion  to  calculate 
strains,  stresses,  and  distortion. 

The  Muraki  program  has  great  flexibility  in  terms  of 
allowing  temperature  dependency  of  material  properties  to  be 
considered  as  well  as  considering  the  effect  of  melting  and 
deposition  of  metal  which  occurs  during  welding.   Weldment 
restraint  and  concentrated  loads  are  also  considered  in  the 
analytical  solution.   In  the  "Manual  on  Axisymmetric  Finite 

Element  Program  for  Analysis  of  Thermal  Stress  and  Metal 

2  8 
Movement  During  Welding",  Muraki   presents  a  brief  synopsis 

of  the  theoretical  background  upon  which  this  analysis  is 

based.   This  synopsis  is  reproduced  here  as  follows: 

Principle  of  Virtual  Work 

The  principle  of  virtual  work  for  the  quasi-static 
problem  is,  in  general, 

///  {<5e}T{a}dV  -  ///  {5u}T{P~}dV  -  //_  {6u}T{F~}dS  =  0    CD 
v  v  1 
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where  superscript  T  stands  for  the  transpose  of  a  vector  or 

•         *  • 

a  matrix.   {e} ,  {a},  and  {u}  are  vectors  of  the  time  rates  of 
strain,  stress  and  displacement,  respectively. 

P  and  F  denote  the  time  rates  of  body  force  over  the 
volume,  V,  and  surface  force  on  the  boundary  surface,  S, , 
respectively. 

The  boundary  condition  on  the  specified  displacement  on 
the  surface,  S?,  is  given  as 

{u}  =  {u}    on  S2  (2) 

It  is  noted  that  the  problem  here  is  formulated  in  the 
small  displacement  theory  based  on  the  rectangular  cartesian 
coordinate  system. 

Constitutive  Equations 

One  of  the  important  features  in  the  analysis  of  thermal 
stress  and  distortion  during  welding  is  the  temperature  de- 
pendency of  material  properties  and  the  yield  criterion. 
Therefore  the  constitutive  equations  are  provided  here  in  a 
general  form. 

The  rate  of  strain  for  the  elasto-plastic  state  is  given 
by: 

ik)    =    {e(e)l  +  {£(p)}  C3) 

where  superscripts  e  and  p  refer  to  elastic  and  plastic  com- 
ponents of  the  strain,  respectively. 
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The  rate  of  plastic  strain  is  assumed  as 

{£(p)}  =  A{g}  (4) 

where  A  is  a  proportional  constant  and  g  is  derivatives  of 

the  yield  function,  f,  with  respect  to  the  stress,  a.. 

(i, j  =  1,2,3) ; 

9f 
9  '  8a±j  (5) 

Using  Hooke ' s  law  and  Equation  (4),  we  have 


{e}  =  ||De|  1    {a}  +  A(g}  (6)* 

where  E|De|  is  the  matrix  form  of  Hooke ' s  law. 

Thermal  strain  and  the  temperature  dependency  of 
material  properties  are  added  to  Equation  (6)  as  follows: 

{e}   =  ElDef1    (a>-  §2-   |De|"1{a} 

(7) 

A  +  e9{5. .}  +  A{g} 

•  9 
where  e   is  the  thermal  strain,  and  6.   is  the  Kronecker 

symbol. 

Equation  (7)  is  rewritten  for  the  rate  of  stress  as, 

{a}  =  E|De|  {e}+  |{a}  -  e9  E|De|{6, .} 

b  13 

-  AE|De|  {g}  (8) 


* 


Note:         stands  for  the  inverse  of  the  matrix 
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If  the  von  Mises1  yield  criterion  is  adopted,  f  in 

Equation  (5)  is  given  by 

f  =  a   -  C  Ce^0  ,  T)  (9) 

where  C  is  related  to  the  strain  hardening  and  given  as  a 
function  of  plastic  strain,  e.?   and  temperature,  T. 

a  is  referred  to  as  the  equivalent  stress  and  given 
by 

a2  =  3/2{a'}T{cT}  (10) 

It  is  noted  that  {cr-*}  is  the  vector  of  deviator  stresses 
which  is  different  from  {a}. 

When  a  state  of  stress  in  a  portion  of  a  body  is  plastic 
and  then  changes  to  a  new  state  of  stress  which  is  still 
plastic,  the  following  relation  is  obtained  from  Equation 
19)  : 

Equation  (11)  with  Equations  (4) ,  (9) ,  and  (10)  leads 
to 

io'}T{'o}    =  -2/3a(-  AH'  +  |f  T)  (12) 


where 


H-  -   -  " 


37(P» 


and  ,-CP»    =  7      2/3{3£lP)}T{3£tP)} 

From  Equations  (8)  and  (12),  A  is  obtained  as  follows: 
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A   = = x  CE{a'}T|Del{e} 

2/3  a  H^  +  E{a'}  |De|{g} 


+  E  {a'}T   {a}  -   e9E{a'}T|De|{6. .} 


+  2/3  a  ||  T)  (13) 


Substituting  A  of  Equation  into  (8)  we  have  a  general 
form  of  the  stress-strain  relationship  as 

{a}  =  (E|De|  -   1/S  E|De|{g}  E{a'}T|De|)  x 
({£}  -  £e{6..})   . 

+  E/E   {|l|    -   1/S      E|De|  {gHa'}T){a} 

°     *e    •  <14) 

-  1/S„  E|De|{g}  2/3a  |£  T 


where 


o   '    '  ^    '    3T 


S   =  2/3  a  H'  +  E{a"}T|De| {g} 


and  |l|  is  called  the  unit  matrix. 

Equation  (14)  covers  the  stress-strain  relationship  in 
the  elastic  as  well  as  plastic  state.   This  can  be  easily 
confirmed  by  eliminating  terms  related  to  the  yield  function 
from  Equation  (14) . 

Yield  Criterion 

As  described  in  Equation  (9)  the  yield  function  used 
here  is  that  proposed  by  von  Mises. 

The  yield  function  is  to  separate  the  region  of  elastic 

deformation  from  that  of  plastic  deformation.   Since  the 
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deformation  is  closely  related  to  the  loading  process,  the 
yield  function  is  also  called  a  loading  surface  which  is 
defined  in  stress  space. 

If  a  part  of  a  material  in  question  is  in  the  region  of 
plastic  deformation,  which  means  that  the  state  of  stress  of 
the  part  is  located  on  the  loading  surface,  one  of  the 
following  states  will  be  realized  at  the  next  instance: 

f  =  0,  f  <  0  (unloading) 

f  =  0,  f"  =  0  (neutral) 

f  a  0,  f  >  0  (loading) 

where 

f   =  *l—a..        2£  t 
3a ij  +  8T 

ij 

If  a  part  of  the  material  is  in  the  region  of  elastic 
deformation,  one  of  the  following  states  will  be  realized 
at  the  next  instance: 

f  <  0    (elastic) 

(16) 
f  =  0    (plastic) 

It  is  worth  looking  into  Equation  (15)  more  closely. 
Using  Equations  (10)  and  (12)  with  Equation  (15) ,  the 
following  relationship  is  obtained: 

f*  =  AH"  (17) 

Equation  (17)  indicates  that  when  non-strain  hardening 
is  assumed  of  a  material,  Equation  (15)  cannot  be  used  to 
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judge  whether  or  not  a  part  of  the  material  recovers  from 
the  plastic  to  the  elastic  state  of  stress.   To  avoid  this 
difficulty,  Equation  (17)  will  be  written  in  the  form  of: 

H 

Therefore  Equation  (15)  is  written  as, 

f  =  0,  A  <  0  (unloading) 

f  =  0,  A  =  0  (neutral)  (19) 

f  =  0,  A  >  0  (loading) 
to  cover  non-strain  hardening  materials. 

Finite  Element  Approach 

Various  types  of  axisymmetric  finite  elements  can  be 

derived  based  on  Equation  (1)  . 

2  8 
In  the  computer  program  described  in  the  Muraki 

Manual,  a  four  node-isoparametric  element  is  used.   Follow- 
ing a  standard  procedure  to  derive  an  element  stiffness 
matrix  and  element  forces,  the  coordinates  and  the  displace- 
ments for  the  elements  are  assumed  by 


x  =  E  h.x. 
.,11 
i=l 


y  =   E   h.y. 
i-1  X  X 


(20) 


u  =   E   h. u.  ,», , 

i=l  x  1  (21) 
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=   Z   h.v. 
i-1  1  X 


(21) 


where  h, ,  h2 ,  h~  ,  and  h4  are  interpolation  functions,  x. 
and  y .  are  nod 
displacements . 


and  y.  are  nodal  coordinates,  and  u.  and  v.  are  nodal 
Ji  li 


A  key  of  the  above  derivation  is  use  of  the  same  in- 
terpolation functions  for  both  the  element  coordinates  and 
the  element  displacements.   Details  of  the  interpolation 
functions  are 


h1  =  htl   -   r)  (1  -  s) 


h2  =  hU  +   r)  (1  -  s) 


(22) 


h3  =  kU   +  r)  (1  +  s) 

h4  =  ha   -  r)  (1  +  s) 

where  r  and  s  are  referred  to  as  natural  coordinates  and  vary 
from  -1  to  +1. 

The  element  strains  in  the  axisymmetric  condition  are 


3x 

3u 

3Y 


u 

Y 
3u 

9y 


A 


(23) 


+ 


3u 

3x 
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Because  u  and  v  in  the  strain  are  a  function  of  x  and 
y,  the  following  relation  is  first  required  between  the  x-y 
and  r-s  systems 


3_ 
9r 


3_ 
8s 


^ 


3x_    3y 


3r     3r 


3x     3y 


9s     3s 


(24) 


Derivatives  of  the  displacements  with  regards  to  r  and 


s  are 


.      x   3h. 
3u  _   „     i 

3r   "   i=19r   Ui 


3u 

3s 


x   3h 


~   =   E 


3u 
3r 


i=l 

x   3h 


3r~  Ui 


(25) 


=   E 


i=l 


3r-  Vi 


_u 
s 


x   h 


v 


i=l 


r   i 


Using  Equations  (24)  and  (25) ,  the  strains  in  Equation 
(23)  are  given  in  terms  of  nodal  displacements,  u,  ,  u« ,  u-., 
and  u.  as  follows: 


{£>  =  |B|{q} 


(26) 


where 
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x    y     z    ' xy 


.rp         ........ 

(q>   =  {u, ,  u  ,  u  ,  u  ,  v  ,  v  ,  v  ,  v  } 


and  {B}  is  the  matrix  which  gives  a  relationship  between 
the  strains  and  the  nodal  displacements.   It  is  noted  that 
{B}  is  a  function  of  r  and  s. 

The  stress-strain  relationship  of  Equation  (14)  is 
rewritten  in  a  simple  form  as: 

{a}  =  |D|{e}  +  {a9}  (27) 

where 

< 

{a}   =  {a  ,  a  ,    a    ,    x      } 
x '   y '   z '   xy 

«  fl 
{D}  is  the  stress-strain  relationship  matrix  and  {cr  }  is  the 

stress  vector  due  to  temperature  change,  and  the  temperature 

dependency  of  material  properties  and  the  yield  function. 

The  element  stiffness  matrix  {k}  is  obtained  using  the 
principle  of  virtual  work  of  Equation  (1)  and  Castigliano ' s 
theorem  as , 

{k}  =   ///{B}T{D}{B}y  d9dxdy  (28) 

27TjV{B}T{D}{B}y|  J|drds 
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where 


J  = 


3x 
3r 

8x 
3s 


3y 
3r 


3y 
3s 


Introducing  a  numerical  integration,  we  have 


{k}  =  2ttZ   {B}7  .{D}{B}  .  .  I  J  I  .  .Y.  .a.   . 

a   a      13     i]1  'i]  i]  id 

-1-  r  J 


(29) 


where  subscripts  i  and  j  indicate  that  Y.  .,  {B} ,  and  |j|  are 
evaluated  at  points  i  and  j.  a. .  is  weight  factor  at  points 
i  and  j  in  the  numerical  integration. 

The  program  in  the  Muraki  manual  uses  the  Gauss-Legendre 
numerical  integration  and  a.   in  that  scheme  can  be  found  in 
text  books  on  numerical  integration. 

The  element  force  vector  {f  }  corresponding  to  the  second 
term  in  the  right  side  of  Equation  (27)  is  obtained  in  a 
similar  manner. 


Tr»e 


{f°}  =  fff{B] x{a  }yd0dxdy 

=   27TjV{B}T{ae}y|j|drds 


(30) 


Again  introducing  the  numerical  integration  we  have 


{fd}  =    2-rrZ       {B}.  .{a°}Y.  .  |j|  .  .a.  . 
•     •         13  3LJ         13    13 


(31) 


In  the  current  version  of  the  computer  program  described 
in  the  Muraki  manual,  only  the  static  load  considered  is  the 
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is  the  concentrated  load  at  nodal  points.   Therefore  there 
is  no  other  element  force. 

The  stiffness  matrix  and  force  vector  for  the  whole 
body  (i.e.  the  whole  system)  are  obtained  by  assembling  the 
element  stiffness  matrices  and  force  vectors. 

The  equilibrium  equations  for  the  whole  system  are 

obtained  as 

Z{k}.{Q}  +  E{f0}.   =   Z{F~}  (32) 

i         i  j    3 

where  i  and  j  stand  for  the  i-th  element  and  the  j-th  node 
respectively. 

{Q}  is  the  vector  of  the  time  rate  of  nodal  displacement 
{F}  is  the  vector  of  the  time  rate  of  concentrated  load. 

Equation  (32)  is  further  written  in  the  incremental 
form  as  follows: 


S{k}i(AQ}  +  E{Af°>i   =   Z(AF}.  (33) 

i  i     x     j     : 


where 

* 

AQ   =   QAt 
Af9  =   f°At 
AF   =   FAt 
At  is  the  time  increment. 

It  is  noted  that  the  time  increment  is  included  in 
Equation  (33)  implicitly. 


93 

Once  the  nodal  displacements  for  the  system,  {AQ},  are 
obtained,  the  element  strains  are  calculated  by  Equation  (26) 
and  then  the  corresponding  stresses  by  Equation  (27) . 

When  there  is  no  further  increment  of  concentrated  and/ 
or  thermal  loads,  the  analysis  is  complete. 

When  there  is  another  increment  of  the  loads,  new 
equilibrium  equations  are  calculated  using  the  nodal  dis- 
placements and  stresses  obtained  in  the  preceding  calculation 
This  procedure  of  incremental  analysis  continues  until  the 
final  increment  is  solved. 


94 

2 .   Geometry  of  Weld  Specimen  and  Finite  Element 
Representation  of  Weldment  used  by  Muraki  Program: 

Finite  element  representation  of  the  weldment  is 
depicted  in  Figure  V-2.   A  comparison  of  Figure  V-2  with 
the  finite  element  representation  of  the  weldment,  Figure  V-l, 
used  in  the  heat  flow  program  will  show  that  an  axis  trans- 
formation as  well  as  a  shift  of  the  origin  were  required  in 
order  to  use  heat  flow  program  output  as  input  for  the  Muraki 
Program. 

For  the  Muraki  program,  the  y-axis  corresponds  to  the 
radial  direction  and  the  x-axis  is  along  the  cylinder  length. 
The  z-axis  corresponds  to  the  direction  of  welding  arc  travel, 
normal  to  the  x-y  plane  as  shown  in  Figure  V-2,  positive 
direction  out  of  the  paper.   The  origin  (x=y=0)  for  the  Muraki 
program  analysis  corresponds  to  a  point  at  the  center  of  the 
cylinder  end.   z=0  corresponds  to  the  time  of  arc  initiation. 

The  Muraki  program  also  required  the  specification  of 
nodal  coordinates  (x,y  values) .   Element  connectivity  was  input 
to  the  Muraki  program  by  identification  of  the  nodes,  in  counter- 
clockwise fashion  beginning  with  the  node  closest  to  the  heat 
source,  associated  with  each  element.   Nodal  coordinates  input 
to  the  Muraki  program  are  listed  in  Table  V-B  and  element  con- 
nectivity is  shown  in  Table  V-C . 

Thermophysical  material  properties  required  by  the 
Muraki  program  were  taken  from  Figures  II-l  through  11-10 
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in  Chapter  II  of  this  study.   Welding  parameters  required  as 

input  data  for  the  first  pass  of  the  girth  weld  were  the  same 
as  those  used  by  the  heat  flow  program.   These  parameters 
were  discussed  and  values  presented  earlier  in  this  chapter. 
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CHAPTER  VI 

Comparison  of  temperatures  and  values  of  strain,  stress, 
and  distortion  during  production  of  the  HY-130  cylinder  girth 
weld  with  those  predicted  by  use  of  the  heat  flow  and  Muraki 
programs  follows  in  the  form  of  tabulated  data  and  in  graph- 
ical presentations.   During  production  of  the  girth  weld,  arc 
travel  began  at  the  seam  weld  opposite  the  cylinder  seam  about 
which  thermocouples  and  strain  gages  were  placed.   Times  at 
which  data  comparisons  are  presented  were  selected  as  repre- 
senting a  thermal  cycle  during  which  temperature,  strains, 
stress,  and  distortion  underwent  greatest  changes. 

Temperatures  determined  experimentally  at  the  eight 
thermocouple  locations  for  the  times  of  interest  are  shown 
in  Table  VI-A.   Distances  from  weld  centerline  are  indicated 
beside  the  column  headings  designated  Tl,  T2,  T3 ,  T4,  T41, 
T3 ' ,  T2 ' ,  and  Tl ' .   These  temperatures  are  presented  graph- 
ically in  Figure  VI-1.   Temperature  distribution  predicted  by 
use  of  the  heat  flow  program  is  presented  by  a  computer  gener- 
ated graph,  Figure  VI-2. 

Qualitative  comparison  of  the  analytical  and  measured 
temperature  distributions  shows  remarkable  coincidence.   The 
slightly  higher  overall  values  of  the  analytically  generated 
values  may  be  explained  by  the  assumed  welding  arc  efficiency 
(narc)  of  0.65. 
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For  the  purpose  of  presenting  comparative  strain,  stress, 

and  distortion  data,  the  following  terms  are  defined: 

STR(*.*)X  =  EXPERIMENTALLY  DETERMINED  STRAIN 

(MICROSTRAIN)  AT  A  POINT  *.*  INCHES 
FROM  THE  WELD  CENTERLINE . 

STR(*.*)A  =  ANALYTICALLY  PREDICTED  STRAIN 

(MICROSTRAIN)  AT  A  POINT  *.*  INCHES 
FROM  THE  WELD  CENTERLINE. 

STS(*.*)X  =  EXPERIMENTALLY  DETERMINED  STRESS  (lbf/.  2) 
AT  A  POINT  *.*  INCHES  FROM  THE  WELD 
CENTERLINE . 

STS(*.*)A  =  ANALYTICALLY  PREDICTED  STRESS  (lbf/.  2) 

7  in 

AT  A  POINT  *.*  INCHES  FROM  THE  WELD 
CENTERLINE. 

DST(*.*)X  =  EXPERIMENTALLY  MEASURED  RADIAL  DISTORTION 
IN  INCHES  AT  A  POINT  *.*  INCHES  FROM  THE 
WELD  CENTERLINE. 

DST(*.*)A  =  ANALYTICALLY  PREDICTED  RADIAL  DISTORTION 
IN  INCHES  AT  A  POINT  *.*  INCHES  FROM  THE 
WELD  CENTERLINE. 

Positive  values  of  strain  and  stress  correspond  to  the 
tensile  condition  and  negative  values  correspond  to  compressive 
strains  and  stresses.  Positive  values  of  radial  distortion 
correspond  to  displacement  in  an  outward  direction  from  the 
cylinder  center.  Tables  and  graphs  on  the  following  pages 
use  the  terms  as  defined  above  to  compare  data.  Strain  and 
Stress  directions  correspond  to  directions  shown  in  Figure  V-2. 
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Tables  VI-B  through  VI-E  present  in  tabular  form 

comparisons  between  strains  measured  during  production  of 
the  HY-130  cylinder  girth  weld  and  those  predicted  by  use 
of  the  Muraki  program.   Tables  VI-F  through  VI-I  present 
the  analytically  predicted  stresses.   Experimentally  deter- 
mined stresses  were  not  tabulated  since  experimentally 
determined  strains  could  not  be  transformed  into  stresses. 
There  are  three  (ignoring  creep)  components  of  strain;  elastic, 
plastic,  and  thermal.   Elastic  and  thermal  constitutive  rela- 
tions are  simple  and  the  stresses  resulting  therefrom  could 
be  calculated  if  these  strains  could  be  separated  into  their 
components.   However,  the  constitutive  relationships  among 
elastic,  plastic,  and  thermal  strains  are  complex,  made  further 
so  by  the  presence  of  a  strain  hardening  material  and  by  the 
dependence  of  all  material  properties  on  temperature.   For 
this  reason,  experimentally  determined  strains  could  not  be 
transformed  into  stresses.   Values  of  radial  distortion  measured 
during  production  of  the  girth  weld  and  predicted  by  the  Muraki 
program  are  shown  in  Table  VI-J.   Comparative  graphs  of  measured 
and  predicted  data  are  shown  in  Figures  VI-3  through  VI-9. 
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LONGITUDINAL  STRAIN 
(X-DIRECTION) 

TRANSVERSE  STRAIN 
(Z-DIRECTION) 

TIME (sec) 

STR(1.0)X 

STR(l.O)  A 

STR(1.0)X 

STR(l.O) A 

125.0 

-44 

0 

200 

0 

127.0 

-130 

0 

235 

0 

128.0 

-153 

0 

277 

0 

129.0 

-177 

10 

318 

18 

129.5 

-190 

-2 

365 

32 

130.0 

-198 

-23 

370 

51 

130.5 

-198 

-127 

405 

132 

131.0 

-185 

-305 

420 

279 

132.0 

-185 

-420 

475 

391 

133.0 

-228 

-350 

607 

387 

134.0 

-390 

-271 

780 

401 

135.0 

-608 

-177 

922 

420 

136.0 

-823 

-68 

967 

446 

137.5 

-933 

127 

902 

481 

140.0 

-818 

474 

527 

531 

145.0 

-464 

1090 

-194 

600 

150.0 

-370 

1530 

-319 

641 

155.0 

-4 

1820 

-184 

668 

160.0 

0 

1930 

-129 

680 

170.0 

49 

2200 

-94 

1 
714   | 

200.0 

176 

2250 

-49 

726 

TABLE  VI-B   EXPERIMENTAL  AND  ANALYTICAL  STRAIN  1.0  IN  FROM 
WELD  CENTERLINE 
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LONGITUDINAL  STRAIN 
(X-DIRECTION) 

TRANSVERSE  STRAIN 
(Z -DIRECT I ON 

TIME (sec) 

STR(1.5)X 

STR(1.5)A 

STR(1.5)X 

STR(1.5)A 

125.0 

-245 

0 

320 

0 

127.0 

-310 

0 

365 

0 

128.0 

-325 

0 

395 

1 

129.0 

-345 

2 

410 

16 

129.5 

-345 

-70 

425 

27 

130.0 

-345 

-21 

450 

41 

130.5 

-340 

-105 

470 

92 

131.0 

-340 

-259 

485 

185 

132.0 

-330 

-358 

505 

260 

133.0 

-330 

-318 

570 

263 

134.0 

-320 

-297 

626 

277 

135.0 

-325 

-288 

681 

293 

136.0 

-330 

-283 

706 

311 

137.5 

-335 

-265 

750 

336 

140.0 

-340 

-205 

749 

373 

145.0 

-330 

11 

664 

433 

150.0 

-350 

276 

560 

477 

155.0 

-400 

527 

496 

513 

160.0 

-390 

651 

427 

530 

170.0 

-330 

1060 

257 

590 

200.0 

-291 

1490 

-246 

666 

TABLE  VI-C   EXPERIMENTAL  AND  ANALYTICAL  STRAIN  1.5  IN  FROM 
WELD  CENTERLINE 
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LONGITUDINAL  STRAIN 
(X-DIRECTION) 

TRANSVERSE  STRAIN 
(Z-DIRECTION) 

TIME (sec) 

STR(2.0)X 

STR(2.0)A 

STR(2.0)X 

STR(2.0)A 

125.0 

-310 

0 

315 

0 

127.0 

-355 

0 

385 

0 

128.0 

-360 

0 

395 

1 

129.0 

-370 

-5 

410 

14 

129.5 

-385 

-16 

435 

22 

130.0 

-400 

-32 

440 

33 

130.5 

-395 

-115 

445 

68 

131.0 

-390 

-264 

455 

132 

132.0 

-380 

-365 

460 

186 

133.0 

-365 

-334 

475 

193 

134.0 

-355 

-321 

490 

207 

135.0 

-320 

-320 

498 

221 

136.0 

-310 

-325 

500 

236 

137.5 

-283 

-331 

502 

257 

140.0 

-225 

-337 

520 

288 

145.0 

-220 

-304 

555 

338 

150.0 

-384 

-216 

606 

378 

155.0 

-509 

-94 

666 

413 

160.0 

-636 

-16 

689 

431 

170.0 

-625 

304 

600 

496 

200.0 

-326 

836 

164 

587 

TABLE  VI-D   EXPERIMENTAL  AND  ANALYTICAL  STRAIN  2.0  IN 
FROM  WELD  CENTERLINE 
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LONGITUDINAL  STRAIN 
(X-DIRECTION) 

TRANSVERSE  STRAIN 
(Z-DIRECTION) 

TIME (sec) 

STR(2.5)X 

STR(2.5)A 

STR(2.5)X 

STR(2.5)A 

125.0 

-175 

0 

295 

0 

127.0 

-220 

0 

310 

0 

128.0 

-250 

0 

320 

1 

129.0 

-265 

-5 

325 

11 

129.5 

-265 

-8 

335 

16 

130.0 

-265 

-13 

355 

24 

130.5 

-265 

-30 

355 

47 

131.0 

-265 

-61 

365 

88 

132.0 

-265 

-86 

365 

124 

133.0 

-260 

-85 

365 

132 

134.0 

-260 

-89 

355 

143 

135.0 

-255 

-93 

355 

154 

136.0 

-245 

-98 

350 

166 

137.5 

-223 

-104 

340 

182 

140.0 

-125 

-114 

340 

205 

145.0 

-115 

-123 

355 

244 

150.0 

-249 

-118 

457 

276 

155.0 

-379 

-99 

552 

305 

160.0 

-501 

-83 

625 

320 

170.0 

-560 

1 

624 

377 

200.0 

-376 

203 

365 

463 

TABLE  VI -E 


EXPERIMENTAL  AND  ANALYTICAL  STRAIN  2 . 5  IN 
FROM  WELD  CENTERLINE 
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LONGITUDINAL  STRESS 
(X-DIRECTION) 

TRANSVERSE  STRESS 
(Z-DIRECTION) 

TIME (sec) 

STS(1.0)X 

STS(1.0)A 

STS  (1.0) X 

STS(1.0)A 

125.0 

0 

0 

127.0 

3 

22 

128.0 

143 

374 

129.0 

436 

593 

129.5 

207 

888 

130.0 

-295 

1236 

130.5 

-2658 

2669 

131.0 

-6614 

5369 

132.0 

-9464 

11431 

133.0 

-8670 

6405 

134.0 

-8263 

5400 

135.0 

"  ' —  "" 

-8235 

3860 

136.0 

-8420 

1908 

137.5 

-8619 

-1588 

140.0 

-8727 

-7665 

145.0 

-7337 

-16930 

150.0 

-4913 

-21859 

155.0 

-2501 

-24190 

160.0 

-1321 

-24732 

170.0 

2508 

-25156 

200.0 

6329 

-21790 

TABLE  VI-F   PREDICTED  X  AND  Z  DIRECTION  STRESSES  1.0  IN 
FROM  WELD  CENTERLINE 
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LONGITUDINAL  STRESS 
(X-DIRECTION) 

TRANSVERSE  STRESS 
(Z-DIRECTION) 

TIME (sec) 

STS (1.5)X 

STS  (1.5) A 

STS(1.5)X 

STS(1.5)A 

125.0 

0 

0 

127.0 

0 

2 

128.0 

7 

29 

129.0 

222 

518 

129.5 

54 

741 

130.0 

-226 

1031 

130.5 

-2336 

1709 

131.0 

-6161 

2950 

132.0 

-8488 

4195 

133.0 

-7210 

' 

4717 

134.0 

-6482 

5298' 

135.0 

-6155 

5751 

136.0 

-6047 

6131 

137.5 

-5846 

6464 

140.0 

-5495 

6320 

145.0 

-4303 

4198 

150.0 

-2819 

1155 

155.0 

-1363 

-1724 

160.0 

-612 

-3085 

170.0 

2010 

-7334 

200.0 

5202 

-10709 

TABLE  VI-G   PREDICTED  X  AND  Z  DIRECTION  STRESSES  1.5  IN 
FROM  WELD  CENTERLINE 
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LONGITUDINAL  STRESS 
(X-DIRECTION) 

TRANSVERSE  STRESS 
(Z-DIRECTION) 

TIME (sec) 

STS(2.0)X 

STS(2.0)A 

STS (2.0)X 

STS(2.0)A 

125.0 

0 

0 

127.0 

0 

1 

128.0 

-6 

20 

129.0 

-13 

382 

129.5 

-293 

487 

130.0 

-716 

608 

130.5 

-3141 

612 

131.0 

-7513 

621 

132.0 

-10343 

974 

133.0 

-9197 

1626 

134.0 

-8632 

2215 

135.0 

-8436 

2661 

136.0 

-8476 

3045 

137.5 

-8508 

3563 

140.0 

-8571 

4237 

145.0 

-8168 

4976 

150.0 

-7258 

4859 

155.0 

-6154 

4135 

160.0 

-5498 

3548 

170.0 

-2939 

887 

200.0 

1115 

-3872 

TABLE  VI-H   PREDICTED  X  AND  Z  DIRECTION  STRESSES  2.0  IN 
FROM  WELD  CENTERLINE 
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LONGITUDINAL  STRESS 
(X-DIRECTION) 

TRANSVERSE  STRESS 
(Z-DIRECTION) 

TIME (sec) 

STS(2.5)X 

STS(2.5)A 

STS (2.5)X 

STS(2.5)A 

125.0 

0 

0 

127.0 

0 

1 

128.0 

-2 

17 

129.0 

-23 

295 

129.5 

-50 

450 

130.0 

-89 

656 

130.5 

-269 

1301 

131.0 

-595 

2465 

132.0 

-827 

3487 

133.0 

-780 

3687 

134.0 

-772 

4001 

135.0 

-782 

4299 

136.0 

-807 

4618 

137.5 

-842 

5036 

140.0 

-889 

5658 

145.0 

-958 

6567 

150.0 

-1025 

7051 

155.0 

-1140 

7139 

160.0 

-1231 

7013 

170.0 

-1689 

5963 

200.0 

-2710 

2510 

TABLE  VI-I   PREDICTED  X  AND  Z  DIRECTION  STRESSES  2,5  IN 
FROM  WELD  CENTERLINE 
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RADIAL  DISTORTION 
(Y-DIRECTION) 

RADIAL  DISTORTION 
(Y-DIRECTION) 

TIME (sec) 

DST(0.5)X 

DST(0.5)A 

DST(1.0)X 

DST(l.O) A 

125.0 

0.0000 

0.0000 

0.0000 

0.0000 

127.0 

0.0000 

0.0000 

0.0000 

0.0000 

128.0 

0.0000 

0.0000 

0.0000 

0.0000 

129.0 

0.0000 

0.0000 

0.0000 

0.0000 

129.5 

0.0000 

0.0002 

0.0000 

0.0000 

130.0 

* 

0.0002 

* 

0.0005 

130.5 

* 

0.0009 

* 

0.0013 

131.0 

* 

0.0026 

* 

0.0028 

132.0 

0.001 

0.0035 

0.001 

0.0038 

133.0 

o.obi 

0.0029 

0.001 

0.0037 

134.0 

0.001 

0.0028 

0.001 

0.0037 

135.0 

0.001 

0.0028 

0.001 

0.0037 

136.0 

0.002 

0.0029 

0.002 

0.0038 

137.5 

0.002 

0.0032 

0.002 

0.0039 

140.0 

0.003 

0.0035 

0.003 

0.0041 

145.0 

0.003 

0.0039 

0.004 

0.0044 

150.0 

0.003 

0.0041 

0.004 

0.0045 

155.0 

0.003 

0.0043 

0.004 

0.0046 

160.0 

0.003 

0.0044 

0.004 

0.0047 

170.0 

0.002 

0.0046 

0.004 

0.0049 

200.0 

0.002 

0.0049 

0.003 

0.0052 

*DEFLECTION  TOO  SMALL  TO  QUANTIFY 
TABLE  VI-J 


EXPERIMENTAL  AND  ANALYTICAL  DISTORTION 
1.0  IN  FROM  WELD  CENTERLINE 
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CHAPTER  VII 

CONCLUSIONS 

This  study  of  residual  stresses  and  distortion  in  high 
strength  steel  weldments  has  shown  that  existing  computer 
programs  require  further  refinement  in  order  to  be  used  as  a 
day-to-day  tool  of  the  welding  engineer. 

The  occurrence  of  residual  stresses  and  distortion  during 
production  of  a  weldment  is  a  complex  phenomena.   Past  efforts 
to  simulate  the  welding  process  involved  clever  simplifying 
assumptions  followed  by  attempts  to  carry  out  experiments  which 
supported  these  assumptions.   The  new  approach  to  welding 
simulation  employs  the  enormous  power  of  computer  techniques 
to  analyze  what  actually  takes  place  during  the  welding  process. 
However,  if  a  large  computational  process  is  used  to  model  the 
complex  welding  process,  the  thermophysical  material  properties 
on  which  the  analysis  is  based  must  be  accurately  known. 

The  axisymmetric  finite  element  program  which  was  developed 
by  Muraki  for  analysis  of  thermal  stress  and  metal  movement 
during  welding  is  promising  in  that  it  showed  good  qualitative 
agreement  with  carefully  generated  experimental  data.   Heat 
flow  and  resulting  temperature  distributions  within  weldments 
can  be  accurately  predicted  by  existing  methods.   Qualitative 
and  quantitative  comparison  between  analytical  and  measured 
temperature  distributions  showed  remarkable  coincidence. 
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CHAPTER  VIII 

RECOMMENDATIONS 

Inasmuch  as  analytical  methods  to  predict  stresses  and 
distortion  resulting  from  thermal  cycles  can  only  be  as  good 
as  the  thermophysical  material  properties  on  which  they  are 
based,  it  is  recommended  that  studies  be  undertaken  to 
accurately  determine  the  temperature  dependence  of  high 
strength  steel  properties  at  alevated  temperatures.   This 
would  enhance  further  efforts  to  simulate  the  welding 
process . 

For  further  study  of  computer  simulation  of  multi-pass 

m 

welding,  it  is  recommended  that  a  micro-processor  be  used 
to  process  welding  parameter  data.   Reduction  of  strains 
requiring  temperature  compensation  is  a  tedious  task, 
especially  if  measurements  are  made  at  several  locations 
and  recorded  on  one  photo  tape.   This  additional  source  of 
error  could  be  reduced  or  eliminated  if  parameters  of 
interest  could  be  calculated  internally  and  displayed 
directly. 

High  temperature  strain  gages  could  significantly  con- 
tribute toward  a  more  meaningful  study.   The  400   F  limitation 
of  gages  used  in  this  study  severely  restricted  their  loca- 
tion. 
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